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Jet-medium interactions are studied via a multi-hadron correlation technique (called "2+1"), 
where a pair of back-to-back hadron triggers with large transverse momentum is used as a proxy 
for a di-jet. This work extends the previous analysis for nearly-symmetric trigger pairs with the 
highest momentum threshold of trigger hadron of 5 GeV/c with the new calorimeter-based triggers 
with energy thresholds of up to 10 GeV and above. The distributions of associated hadrons are 
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studied in terms of correlation shapes and per-trigger yields on each trigger side. In contrast with 
di-hadron correlation results with single triggers, the associated hadron distributions for back-to- 
back triggers from central Au+Au data at ■ v /ijviv=200 GeV show no strong modifications compared 
to d+Au data at the same energy. An imbalance in the total transverse momentum between hadrons 
attributed to the near-side and away-side of jet-like peaks is observed. The relative imbalance in the 
Au+Au measurement with respect to d+Au reference is found to increase with the asymmetry of the 
trigger pair, consistent with expectation from medium-induced energy loss effects. In addition, this 
relative total transverse momentum imbalance is found to decrease for softer associated hadrons. 
Such evolution indicates the energy missing at higher associated momenta is converted into softer 
hadrons. 

PACS numbers: 25.75.-q 



I. INTRODUCTION 



Angular di-hadron correlations with respect to a sin- 
gle charged or neutral high-px trigger at the center-of- 
mass energy per nucleon pair v /s]v r /v= 200 GeV have been 
found to differ significantly between heavy ion events and 
more vacuum-like pp or d+Au collisions. On the away- 
side (Aip ~ 7r) of the trigger hadron, broader correla- 
tion distributions and softer transverse momentum (pt) 
spectra of associated hadrons have been reported for cen- 
tral Au+Au events [H, . In some associated hadron pt 
ranges (~2 GeV/c) the modified away-side no longer re- 
sembles the jet-like peak but shows a concave shape near 
7T. Novel features have also been discovered in Au+Au 
data on the near-side (small relative azimuth) of the trig- 
ger hadron: a long-range longitudinal plateau in rela- 
tive pseudorapidity (A?y), called the "ridge" 0. Multiple 
theoretical models have been proposed to simultaneously 
explain these structures, including in-medium parton en- 
ergy loss [|[ and initial state medium fluctuations lead- 
ing to higher-order flow components such as the triangu- 
lar flow (v 3 ) @. However, experimental measurements 
based on 2-particlc correlations with respect to high-p^ 
trigger introduce surface bias for the initial hard scat- 
tering. Thus the partons and subsequently formed jets 
on the near- and away-sides can be affected by differ- 
ent underlying physics. On the other hand, di-hadron 
correlations between two high-pT particles exhibit jet- 
like peaks in both near- and away-sides [H, 0] with little 
shape modification from d+Au to central Au+Au, but 
a strong suppression on the away-side amplitude. This 
observation may be interpreted in two jet-medium inter- 
action scenarios. One possibility is in-medium parton 
energy loss followed by in-vacuum fragmentation, which 
naturally explains the observed high-px inclusive hadron 
suppression, Raa, defined as a ratio of spectra measured 
in Au+Au collisions with respect to binary-scaled pp ref- 
erence @, Q . Such scenario can also explain little to no 
modifications of near-side peaks in the two-particle cor- 
relations (Iaa) [1,0; however, lack of broadening on the 
away-side of such correlation functions poses a challenge 
for this model Q. Alternatively, there could be a finite 
probability for both partons to escape the medium with- 
out interactions. In the model this could be simplified as 
the "core/corona" model, where the dynamics of jets and 



di-jets in the medium are the same as in vacuum, unless 
they traverse the "core" of the medium where they are 
fully absorbed. Varying the surface bias experimentally 
could allow differentiating between the models. 

The analysis presented in this paper uses a 3-particle 
("2+1") correlation technique for the data from the 
Solenoidal Tracker at RHIC (STAR) [nj. A pair of 
back-to-back high-py trigger particles is used as prox- 
ies for di-jet axis, and angular correlation of lower-pT 
charged hadrons with respect to back-to-back trigger pair 
is considered. This technique was first introduced in 
the previous STAR publication [H| . In the early work, 
the primary trigger (Tl) and its back-to-back secondary 
trigger partner (T2) had similar kinematic thresholds of 
(p^ 1 > 5 GeV/c and plj? > 4 GeV/c). Following previ- 
ous works, the near-side is defined as rj — (f> region close 
to trigger Tl, and away-side is the r\ — 4> space close to 
the away-side trigger T2. It has been found, that such 
di-jet-like correlations from central Au+Au collisions are 
similar in both shape and magnitude to those observed in 
d+Au events at the same incident energy on both near- 
and away-sides. This suggested a strong surface (tangen- 
tial) bias of selected di-jets |11|. 

In this paper the analysis is further extended to higher 
primary trigger thresholds. We attempt to use the asym- 
metry in the energy of the two back-to-back triggers as 
a tool for changing the surface bias. The expectation is 
that such asymmetry can partially arise from a longer 
path length that the away-side parton must travel in the 
medium, thus producing differences in the balance of fi- 
nal jet energies. To control the degree of the surface bias 
we vary the relative balance between the energies of the 
primary trigger Tl (near-side) and its back-to-back trig- 
ger partner T2 on away-side. For the most asymmetric 
trigger pair selection in this work the primary trigger has 
more than twice the energy of its back-to-back partner. 
The correlation functions and the spectra of associated 
charged particles around each trigger within A0 < 0.5 
and |At/| < 0.5 are measured and compared between 
Au+Au and d+Au reference data. 

The total transverse momentum of each side of a di-jet 
is then calculated by summing the px for all associated 
charged hadrons plus the trigger Et or px- The di-jet 
energy imbalance, A(EEt), is calculated as difference be- 
tween the total transverse momentum between the same- 
side and away-side jet-like peaks. The absolute value of 
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the imbalance could be affected by the kinematic selec- 
tion of the trigger pair unrelated to the jet- medium inter- 
actions, for example, by the kx effect. To provide quan- 
titative assessment of jet-medium interaction effects and 
allow discriminating of the theoretical models, the rela- 
tive di-jet energy imbalance between Au+Au and d+Au 
data is more informative. 



II. DATA SETS 

In this work the d+Au and Au+Au collisions at 
Y / J/vw=200 GeV recorded by the STAR collaboration in 
Runs 2007 and 2008 are analyzed, extending analysis of 
the previous, smaller, samples from Runs 2003 and 2004. 



Earlier measurement has been previously reported in [11 



we will use these results here for comparison. The details 
of new datasets and analysis selections are discussed be- 
low. For a more uniform detector acceptance only events 
with a primary vertex position Vz within 25 cm of the 
center of the STAR Time Projection Chamber (TPC) 
along the longitudinal beam direction were used in the 
analysis. Run 2007 provided 74M Au+Au minimum-bias 
(MB) events. These MB events were selected by requiring 
at least one hit in the Vertex Position Detectors (VPD) 
located on each side of the TPC 4.5 m away from the 
nominal interaction point. A subset of such events with 
additional requirement of at least one high-energy tower 
(with transverse energy Et > 5.75 GeV) in the Barrel 
Electromagnetic Calorimeter (BEMC) is referred as the 
"high-tower trigger data" . Run 2008 collisions provided 
6M d+Au high-tower triggers (with Et > 4.3 GeV) and 
46M VPD minimum-bias triggered events. To maximize 
any medium effects, this work focuses on the analysis of 
0-20% most central Au+Au events. 

The BEMC clusters are each built from a set of closely 
neighboring calorimeter towers and showcr-max detector 
strips with appropriate shape and quality cuts. BEMC 
cluster with the highest transverse energy of at least 
8 GeV was selected as the primary trigger (Tl). The trig- 
gered events were grouped into two different bins with 
E T £ [8, 10] GeV and E T G [10, 15] GeV. A highest- 
Pt charge particle in the back-to-back azimuthal region 



(\0 



Tl 



xT2 



7r| < 0.2) is selected as the back-to- 



back partner trigger (T2) with kinematic requirement 
of 4 < pP < 10 GeV/c. This width of 0.2 is an ap- 
proximate width of the away-side peak in the high-pr 
trigger-trigger correlation Q, ensuring a good balance 
between the di-jet purity and signal-to-background ratio. 
For proper primary trigger designation it is also required 
that the momentum of the secondary trigger is less than 
that of the primary one. The new data sets, especially 
benefiting from implementation of triggering capabilities 
of the STAR BEMC detector, allowed for significant im- 
provement of the kinematic reach of this study, and thus 
larger asymmetry span in the back-to-back trigger part- 
ner selection. The primary trigger Tl threshold has been 
moved up by a factor of two while still maintaining a 



good energy resolution. The correlation is constructed 
with all other charged hadrons in an event. These associ- 
ated particles are required to have transverse momentum 
1.0 < pjp soc < 10 GeV/c, and are divided into groups 
with different lower thresholds, [1,10], [1.5,10], [2,10], 
and [2.5, 10] GeV/c. The kinematic selection for associ- 
ated hadrons is intended to cover the region where the 
broadening of the away-side correlation and the near-side 
ridge were previously reported in the 2-particle correla- 
tion measurements 0, H| . 

The associated charged particles reconstructed by 
TPC are required to pass a set of track quality cuts. Mid- 
rapidity 1 77 1 < 1.0 tracks selected are required to have at 
least 20 fit points for good momentum resolution, and a 
distance of closest approach (DCA) to the primary ver- 
tex of less than 1 cm. A fiducial cut of |r/| < 0.75 is also 
required for the center of each EMC cluster. The high- 
Pt charged particle contamination in the BEMC trig- 
ger sample is removed by the following charged veto cut: 
a BEMC-cluster trigger is rejected if a charged parti- 
cle with pt > 1.5 GeV/c is projected to a point within 
\A(f>\ < 0.015 and [A^l < 0.015 of the center of the clus- 
ter. 



III. ANALYSIS METHOD 

A. Building the Correlation Functions 

The correlation functions studied in this work are de- 
fined as 



d 2 N 



1 



d 2 N I£ 



dArjdAip N tl i g e pa i r dAr)dA<f) 



-) 



x zyam 



dA?]dA(j) 



(1) 



where Ntng is the number of trigger pairs, 
d 2 N lav/ / dAr/ dAcf> is the associated hadron distribu- 
tion relative to each trigger in correlated (Tl, T2) pair, 



and 



pair 



is a correction factor for single-track efficiency 



and pair acceptance effects. The d 2 N-Q g /dArjdA(j> rep- 
resents a background term, originating predominantly 
from randomly associated pairs and correlations due to 
anisotropic flow. Background scaling factor, a zyam , is 
described in the text below. The tracking reconstruction 
efficiency is derived by embedding single Monte- Carlo 
tracks into the real events and reconstructing the 
combined event. To account for the 77, 0, pt and 
multiplicity dependence of the single track reconstruc- 
tion efficiency, the correction factor is calculated for 
each track individually and applied on a track-by-track 
basis. The pair acceptance correction is derived by 
the mixed-event technique. For the event mixing all 
accepted events with primary vertices \Vz\ < 25 cm are 
divided into ten 5 cm- wide bins. The events are also 
grouped into 3 centrality sub-classes, corresponding to 
0-5%, 5-10%, and 10-20%. The triggers arc then only 
mixed with associated particles from minimum bias 
events of the same Vz and centrality bin to better reflect 
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FIG. 1: (Color online) (a) A sample of two-dimensional raw correlation function from 200 GeV 0-20% central Au+Au data before 
efficiency and acceptance correction. Primary trigger Ej- 1 € [8,10] GeV, secondary trigger -pTj? £ [4,10] GeV/c, associated 
hadrons 1.5 < p?j? s ° c < 10 GeV/c. (b) Corresponding two-dimensional mixed-event pair-acceptance correction, (c) Same 
correlation function corrected for efficiency and pair acceptance effects. 



acceptance of real events and avoid potential trigger 
biases. The two-dimensional (Arj-A<j>) mixed-event 
correlation function is scaled such that the highest 
Arj ~ bin is normalized to unity (see e.g. Fig. [T]b). 
The raw correlation function (iV raw ) is then divided 
by this normalized mixed-event distribution, shown in 

Fig. in 



B. Background Subtraction 



The background term 



d z N Bg 



in EqJT] for each cor- 



dAr]dA<p 

relation function originates predominantly from random 
combinatorics and correlations induced by collective flow. 
The shape of this background is described by multiple 
components of a Fourier decomposition, with main con- 
tribution in the kinematic region of this analysis from the 
second-order Fourier component usually associated with 
elliptic flow (i>2). The multiplicity and pt dependence of 
i>2 for triggers and associated particles are obtained from 
existing STAR measurements [l2j]. To reduce the effect 
of fluctuations in our v 2 estimate we average the results 
of Event Plane and Four-particle Cumulant methods [l3[ . 
For pt > 4 GeV/c the elliptic flow magnitude is assumed 
to be constant and at the level reported in the high-p^ 
pion flow measurement [3] ■ This assumption is justified 
since the majority of charged hadrons at these momenta 
are charged pions [1, Q, and high-£r clusters used in 
this work are mainly produced by high-py tt° (HI- Due 
to the back-to-back requirement of trigger pair selection, 
the distribution (apart from efficiency and acceptance ef- 
fects) is modulated by f(A(f>) = l+v 2 +1 cos(2A(/>), where 
i>2 +1 is resulting flow modulation for three- particle cor- 
relation [H, [I?} given by 



,,2+1 



Q„,T1 or T2„, assoc 
£U 2 U 2 



2v 



T2 or Tl„, assoc sin(2c 



la 



1 + 2^X2^ 



(2) 



Here a=0.2 is the half-width of the back-to-back trigger 
cone. The overall background level a zyam is estimated 



with the Zero- Yield at Minimum method [18l - l20l |. Each 
2-D correlation function within |A7j| < 1.0 is first pro- 
jected on relative azimuth to optimize signal-to-noise ra- 
tio and avoid fluctuations at the edges of the TPC ac- 
ceptance. Then the zero-yield region for the "2+1" cor- 
relation is chosen to be consistently at least 1.3 radians 
away from both jet-like peaks at Acf> = or tt . We note, 
that this is more than 3c of jet-like peak widths if fit 
by a Gaussian. A double Gaussian plus a ^-modulated 
background fit is also used for the background level es- 
timate to evaluate systematic uncertainty of the ZYAM 
method. The transverse momentum spectra for associ- 
ated hadrons in the jet peaks is obtained in a similar 
manner from the pT-weighted correlations selecting the 
hadrons within 0.5 radians in relative azimuth and 0.5 in 
relative pscudorapidity of the respective trigger direction. 

An additional background term is related to the ran- 
domly associated triggers in the initial selection of the 
trigger pairs (T^ |. The signal-to- noise ratio is measured 
from the trigger-trigger correlation to estimate the rela- 
tive contribution per trigger pair for such random asso- 
ciations. The correlation contribution due to this back- 
ground is constructed from two independent 2-particle 
correlations for Tl and T2 trigger selections separately 
as in llj. The ZYAM method is also applied to these 
2-particle correlations, and the zero-yield region is set 
to be [0.8, 7r — 1.9] relative to the trigger, also to avoid 
both jet-like peaks on near- and away-sides. An example 
trigger-trigger correlation to illustrate signal to noise ra- 
tio in trigger pairs is presented in Fig. [2] a. The dashed 
line shown in the figure illustrates the ZYAM level. The 
number of true di-jet triggers and random pairs is ob- 
tained from counting the entries above and below the 
ZYAM, respectively. The 2-particle correlation used as 
an estimate of correlated background contribution corre- 
sponding to kinematic selection of first and second trig- 
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gers are shown in Fig. [2] b and c, respectively. We note 
the known ridge- like structures evident in these 2-particle 
correlations as expected. 

A sample set of fully corrected "2+1" correlation func- 
tions is shown in Fig. [31 with the background induced 
by uncorrelated trigger pairs and the ^-modulation re- 
moved. In the following sections we will present all the 
correlation functions measured for the small relative an- 
gles only on both trigger sides to focus on the jet-like 
peaks. 

An additional test was conducted to estimate possible 
effects from higher-order flow terms, such as v^, on the 
correlations. This was done by including higher-order 
flow terms into the fit of final correlation functions such 
as in Fig. [3l We found the effects of v% to be negligible in 
the kinetic regime of this analysis. This result is reason- 
able as higher-order flow terms are generally explained by 
the "soft" and "non-jet" sources; thus the jet-dominant 
source of our trigger pairs could naturally lead to little 
or no v 3 contribution. Independent measurements plj 
also confirm that the contributions from higher-order i; 3 
and U4 terms are much smaller than v 2 in the trigger pt 
region of this paper. Therefore no corrections for flow 
terms other than v 2 are included in this analysis. 



IV. RESULTS 
A. Correlation Functions and Spectra 

In this section we present results obtained from fully 
corrected correlations for each pr bin studied for both 
d+Au and Au+Au datasets. The 2-D correlations them- 
selves are shown in Figs. |4][7] for primary triggers with 
transverse energy 8 < E^ 1 < 10 GcV and 10 < E^ 1 < 
15 GeV. The secondary back-to-back trigger for all cor- 
relations shown is chosen to have 4 < p*P < 10 GeV/c. 
Only statistical errors are included in each of these plots. 

In Figs. [S] and |H] the A</> and Ar/ projections of 2-D 
correlation plots for the asymmetric back-to-back triggers 
studied are shown. The associated hadron transverse mo- 
mentum selection for projections shown is [1,10] GeV/c, 
the lowest associated pr bin studied. Also shown are the 
transverse momentum spectra for the associated charged 
particles within a 0.5x0.5 (A77 x A</>) area of each trig- 
ger, the region containing the dominant part of the jet- 
like peak. Only statistical errors are included in each of 
the plots; the systematic errors are strongly correlated 
between the near- /away-side signals of each collision sys- 
tem and are discussed in detail in the text. We find 
that for each Tl bin studied, the correlation functions 
and associated particle spectra from Au+Au events are 
similar to those observed in the d+Au data in spite of 
increased asymmetry between Tl and T2 triggers. The 
near-side yields of both the d+Au and the Au+Au data 
are smaller than their respective away-side yields, and 
the near-side Au+Au yield slightly drops when Et of Tl 
increases from [8,10] to [10,15] GeV. Both observations 



are possibly due to the increased direct 7 contamination. 
This effect will be further discussed later. 



B. Systematic Errors 

Several sources of systematic uncertainties have been 
evaluated as outlined below. The single track recon- 
struction efficiency for each centrality and px bin was 
derived as a 2-D function of rj and <j). The estimated 
uncertainty for this correction is 5% in each centrality- 
PT-rj-(f> bin. This uncertainty is reduced to 1% when the 
correlated yields of the near-side are compared to those 
of the away-side from the same data sample since the 
multiplicity-dependent effects cancel. The finite statisti- 
cal uncertainty on the mixed-event correlation leads to 
a systematic uncertainty in the pair-acceptance correc- 
tion. This error is estimated to be less than 5% for each 
centrality bin. The uncertainty due to anisotropic flow 
contribution is estimated from the difference in the v 2 re- 
sults from the Event Plane and Four-particle Cumulant 
methods. While the fluctuations of v 2 itself can be large, 
the final effect is relatively small due to the evident jet- 
like peak shapes on both trigger sides. This uncertainty is 
found to be 5% in the central Au+Au. This uncertainty 
is largely correlated between the near- and away-sides. It 
is not applicable to d+Au events where no v 2 modulation 
is included. Since higher-order Fourier terms, such as W3 
and U4, have little effect on this analysis, no systematic 
uncertainties are assigned from these sources. 

The systematic uncertainty due to the ZYAM normal- 
ization of the background level was estimated by vary- 
ing the Atp range from which the minimum used in the 
ZYAM method was derived. Other background assump- 
tions, such as a double-Gaussian plus ^-modulated back- 
ground are also used to estimate the scale of this un- 
certainty. The corresponding systematic uncertainty is 
found to be less than 5% in this analysis. This uncer- 
tainty is strongly correlated between near- and away- 
sides and cancels in relative side-to-side comparisons. 

The systematic uncertainty due to the correlated back- 
ground subtraction from the di-jet sample is determined 
to be less than 3% for both d+Au and Au+Au events, 
and is particularly small in the high-£/r trigger ranges 
where signal-to-background ratio is high. This error is 
also correlated between the near-/away-sides, and is esti- 
mated by varying the background normalization for the 
trigger-trigger correlation in a manner similar to that 
used for the "2+1" correlation. Since 2-particle cor- 
relations were used for the correction of random com- 
binatorics of trigger pairs, the uncertainty in such di- 
hadron distributions arises mainly from the aforemen- 
tioned sources: uncertainty in v 2 and ZYAM normal- 
ization. These uncertainties were evaluated in a similar 
manner as above and the results are estimated to vary 
from less than 1% in d+Au events up to about 5% in 
central Au+Au data. 

Despite the charged-veto cut, there are still possible 
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FIG. 2: (Color online) (a) Azimuthal trigger-trigger correlation from 0-20% central 200 GeV Au+Au collisions. Primary trigger 
E} 1 £ [8,10] GeV, secondary trig ger py 2 £ [4,10] GeV/c. Dashed line illustrates the ZYAM level. The contribution from 
randomly associated pairs for this analysis is visualized by the green area; from true di-jets - by the red one. (b) Di-hadron 
correlation for charged hadrons with 1.5 < p^ soc < 10 GeV/c with respect to triggers matching Tl selection, (c) Di-hadron 
correlation for charged hadrons with 1.5 < p5? soc < 10 GeV/c with respect to triggers matching T2 selection. 
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FIG. 3: (Color online) Fully corrected "2+1" correlation function from 0-20% central Au+Au events with respect to the primary 



trigger (a) and secondary trigger (b). 
hadrons 1.5 < p?" < 10 GeV/c. 



Primary trigger £ [8,10] GeV, secondary trigger £ [4,10] GeV/c, associated 



contributions to the BEMC cluster trigger energy from 
\ow-px charged particles. This contribution is estimated 
via averaging the energy of those BEMC towers which 
are far from each trigger and pass the hot map cuts. The 
charged track signal contamination to the high energy 
tower cluster was estimated to be 2.5% in central Au+Au 
and less than 1% in d+Au. 

The overall systematic uncertainty in each trigger- 
associated combination bin is less than 15% after all the 
contributions mentioned above are summed in quadra- 
ture. This is the level of uncertainties important for ab- 
solute measurements of the near-/away-side correlation 
functions or spectra in each collision system. The sys- 
tematic errors are dominated by sources strongly corre- 
lated between near-/away-sides, and will therefore mostly 
cancel when near- and away-sides are compared. Thus, 
relative measures, such as energy imbalance, have higher 
sensitivity to physics effects. 



V. DISCUSSION 

Our previous results [ll| have shown similarity of the 
jet-like peaks not only on the near-side, but also on the 
away-side of the primary high-py trigger between central 
Au+Au and d+Au data for the hadrons in the kinematic 
range of pr > 1.5 GeV/c associated with the back-to- 
back high-pT trigger pair. In addition, no evidence of 
the "dip" or "ridge" measured in di-hadron correlations 
with respect to a single trigger of similar kinematic selec 
tion is present in the "2+1" correlations reported in 11 



This similarity in the correlation shapes was further sup- 
ported by the similarities in the associated hadron px dis- 
tributions. Significant softening of the away-side spectra 
observed in di-hadron correlations [l|, one of the known 
indications of energy deposition into the medium, was 
not evident in these "2+1" data. A simple accommoda- 
tion for the observed differences between di-hadron and 
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FIG. 4: (Color online) Near-side associated hadron distributions about primary triggers from 200 GeV d+Au data. The top row 
shows the correlations for the primary trigger Et selection of [8,10] GeV, the bottom row - for [10,15] GeV. For all correlations 
the away-side trigger pr is in [4,10] GeV /c. From left to right the associated pr ranges are [1,10] GeV/c, [1.5,10] GeV/c, 
[2,10] GeV/c, and [2.5,10] GeV/c. 




FIG. 5: (Color online) Away-side associated hadron distributions about secondary triggers from 200 GeV d+Au data. The 
top row shows the correlations for the near-side trigger Et selection of [8,10] GeV, the bottom row - for [10,15] GeV. For 
all correlations the away-side trigger pr is in [4,10] GeV/c. From left to right the associated pr ranges are [1,10] GeV/c, 
[1.5,10] GeV/c, [2,10] GeV/c, and [2.5,10] GeV/c. 



"2+1" correlation results could be provided by the tan- 
gential di-jet emission scenario. In the presence of very 
strong energy loss in the core of the medium, only back- 
to-back parton pairs produced close to the surface and 
traversing the minimal amount of the medium could be 



recovered in the analysis. The Au+Au correlation func- 
tions would then naturally become similar to those of 
d+Au. The "2+1" results from symmetric triggers do 
not, however, exclude path-length dependent energy sce- 
nario, where finite in-medium energy loss is followed by 
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FIG. 6: (Color online) Near-side associated hadron distributions about primary triggers from 200 GeV 0-20% central Au+Au 
data. The top row shows the correlations for the near-side trigger Et selection of [8,10] GeV, the bottom row - for [10,15] GeV. 
For all correlations the away-side trigger pr is in [4,10] GeV/c. From left to right the associated pr ranges are [1,10] GeV /c, 
[1.5,10] GeV/c, [2,10] GeV/c, and [2.5,10] GeV/c. 




FIG. 7: (Color online) Away-side associated hadron distributions about secondary triggers from 200 GeV 0-20% central Au+Au 
data. The top row shows the correlations for the away-side trigger Et selection of [8,10] GeV, the bottom row - for [10,15] GeV. 
For all correlations the away-side trigger pr is in [4,10] GeV/c. From left to right the associated pr ranges are [1,10] GeV/c, 
[1.5,10] GeV/c, [2,10] GeV/c, and [2.5,10] GeV/c. 



in-vacuum fragmentation, for all di-jets. By increasing 
the asymmetry in transverse momentum between the two 
back-to-back trigger hadrons one may control the degree 
of the surface bias, and thus the in-medium path-length 
traversed by the parton. The path-length effects on en- 



ergy loss then could be studied through the imbalance 
induced by the jet-medium interaction on the final ener- 
gies of each side of di-jet. 

In this work the asymmetry between the back-to-back 
triggers is exceeding a factor of two, however, the corre- 
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FIG. 8: (Color online) Projections of 2-D correlation functions on A(/> (a) (with |At?| < 1.0) and At] (6) (with \A</>\ < 0.7) 
for the hadrons associated with their respective triggers (Tl for near-side, T2 for away-side) are shown for d+Au (circles) 
and central 0-20% Au+Au (squares). Errors shown are statistical. The kinematic selection is as follows: 8 < E^ 1 < 10 GeV, 
4 < p£ 2 < 10 GeV/c, 1.0 < py soc < 10 GeV/c. c) Transverse momentum spectral distributions per trigger pair for the near- 
and away-side hadrons associated with di-jet triggers {\Acj>\ < 0.5, |Ar;| < 0.5). The inclusive charged hadron distribution from 
the central 0-10% Au+Au data is shown for comparison. 



Asymmetric Trigger II: T1 s [10,15] GeV, 12 e [4,10] GeV/c 
■ — Au+Au, Same-side — • — d+Au, Same-side — a — Au+Au, Away-side — © — d+Au, Away-side 




FIG. 9: (Color online) Projections of 2-D correlation functions on Acf> (a) (with \Arj\ < 1.0) and Arj (b) (with \A<f>\ < 0.7) for 
the hadrons associated with their respective triggers (Tl for near-side, T2 for away-side) are shown for d+Au (circles) and 
central 0-20% Au+Au (squares) data. Errors shown are statistical. The kinematic selection is as follows: 10 < E^ 1 < 15 GeV, 
4 < p^ 2 < 10 GeV/c, 1.0 < py soc < 10 GeV/c. c) Transverse momentum spectral distributions per trigger pair for the near- 
and away-side hadrons associated with di-jet triggers (\A(f>\ < 0.5, IA77I < 0.5). The inclusive charged hadron distribution from 
the central 0-10% Au+Au data sample is shown for comparison. 



lations still show similar jet-like peak shapes and magni- 
tudes from d+Au to central Au+Au collisions for both 
the near- and away-sides. No evident ridge or dip struc- 
ture is observed at cither near- /away-side in the Au+Au 
data. We point out that statistical limitations prevent 
the complete exclusion of a ridge which could have 
a very small magnitude in the high-p^ trigger regime 
studied here. The absence of strong modifications in the 
Au+Au correlations relative to d+Au signifies the inten- 
sity of jet-medium interactions, as the recovered di-jcts 
still appear to have significant surface bias even for the 
large trigger asymmetry. We note an increase in the cor- 
relation magnitude on away-side with respect to near- 
side in both Au+Au and d+Au correlations for our most 
asymmetric trigger pair selection. This increase by it- 



self, as observed in both systems, is likely not related 
(or not in full) to the parton/jet energy loss. Since the 
trigger pairs pre-selected some di-jet asymmetry, the to- 
tal yields and Spy of the associated particles from the 
higher-cnergy-trigger side must be lower due to energy 
conservation without any energy loss. Hence the relative 
change between Au+Au and the reference d+Au data is 
a more relevant measure. Additionally, since the major- 
ity of the systematic errors are strongly correlated in this 
analysis between the near- and away-sides, the measure- 
ments of relative differences also have an advantage from 
experimental point. The spectra of associated charged 
particles for each trigger side from central Au+Au events 
are compared to the corresponding distributions from 
d+Au data in Fig.[5]and Fig.O The ratios of Au+Au to 
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d+Au spectra are close to unity and have no prominent 
Pt dependence, indicating that the softening or suppres- 
sion, if present, is not nearly as strong as observed in the 
single particle Raa- 

An additional factor to consider in the comparisons is 
a possible direct-7 contamination of the primary triggers 
measured by BEMC clusters, since the near-side yields of 
direct-7 triggers are expected to be close to zero 0, |23| , 
potentially reducing the per-trigger measured yields on 
the near-side. While the primary triggers from BEMC 
clusters of 8 < Ej, 1 < 15 GeV are dominated by pho- 
ton pairs from ir° decays, grouped into the same BEMC 
cluster due to large tower size, the direct 7 contamination 
may not be negligible over the Et range of this study [HI . 
We expect, however, that such direct-7 contamination 
would reduce the energy imbalance measured if affecting 
it at all: since no in-medium energy loss is expected for 
direct-7, the away-side jet corresponding to such trig- 
ger would be relatively higher than that of 7r° trigger 
due to the away-side (reverse) surface bias. Meanwhile, 
the consistent magnitudes of the near-side Au+Au and 
d+Au correlation functions indicate that the contribu- 
tion of direct-7 triggers is relatively small. In short, the 
"2+1" correlations can be used as tool for more differ- 
ential measurements of path-length effects on in-medium 
energy loss for hard scattered partons. 

The difference between the near- and away-side jet en- 
ergy is used to quantify the medium effects. The jet 
energy is estimated by summing the Et (pt) of trigger 
and associated charged particles within the 0.5x0.5 (Ar] 
x Atfi) area of each trigger in the kinematic range used 
in this analysis. The imbalance between the near- and 
away-sides is measured by the difference of the corre- 
sponding energy sums: 

A(££ T ) = (E^ 1 + £ p ^f° c ' near ) _ ( p T2 + Ep -oc,awa y) 

This energy imbalance, A(EEt), has smaller systematic 
uncertainty than energy sums for the correlated peaks 
themselves, as most sources of systematic errors cancel 
out as discussed above. A non-zero value of A(Y,Et) 
could result from QGP medium effects and/or the known 
/cT-effects characteristic for back-to-back partons [24[ • To 
disentangle these two contributions the measured energy 
imbalance for di-jets from central Au+Au events is com- 
pared with the corresponding value obtained from d+Au 
data with the same trigger/associated particle px cuts. 
The energy imbalance for symmetric trigger pairs (llj 
was found to be 1.59 ± 0.19 GeV/c in central Au+Au 
collisions, similar to the value of 1.65 ± 0.39 GcV/c re- 
ported for the minimum bias d+Au data. This value is 
also close to the initial state kj_ effects of ~1.6 GeV/c es- 
timated in di-jet correlation [J] and disfavors additional 
partonic energy loss into the medium for these Au+Au 
data. The detailed results for our new measurement for 
the asymmetric trigger pairs are summarized in Tables 
U and |TTJ In each table, only the statistical uncertainties 
are listed for each variable. The statistical uncertainties 
on mean values of trigger Et or pt are negligible. 



In Fig. 1101 the relative energy imbalance between cen- 
tral Au+Au and d+Au data (e.g. A(££r) Au+Au - 
A(££r) d+Au ) is shown for both the symmetric (solid 
symbol) and asymmetric trigger cases (open symbols). 
The errors shown include both statistical and systematic 
uncertainties discussed in the previous section. 




Lower threshold of p"*° c (GeV/c) 

FIG. 10: (Color online) The relative di-jet energy imbalance 
estimate A(E£ T ) Au+Au - A{EE T ) d+Au . The open circles 
show results for primary triggers in [8,10] GeV, the open 
squares - [10,15] GeV. The solid data point is the result for 
symmetric trigger pairs from previous STAR paper [Tl|] . The 
x-axis shows lower threshold of the transverse momentum 
selection for associated hadrons. 

The observed relative energy imbalance 
A(E£; T ) Au + Au - A(££r) d+Au between Au+Au 
and d+Au for the asymmetric trigger pairs studied is 
consistent with non-zero values. The relative imbalance 
is consistently higher for all associated hadron selections 
for more asymmetric triggers, as expected if larger 
asymmetries are related to longer paths traversed and 
higher energy loss. The measured values of relative 
energy imbalance are significantly smaller than the 
theoretical prediction for in-medium energy deposition 
based on path-dependent energy loss In this model 
approximately 3 GeV is expected for our asymmetric 
case. The "2+1" results for both symmetric and 
asymmetric trigger pairs seem to favor much stronger 
surface bias than expected in theory and thus point to 
stronger path-length dependence of energy loss. Also, 
as previously discussed, any direct-7 contamination to 
our primary trigger sample would reduce the observed 
energy imbalance, leading to even greater discrepancy 
between the data and the theory. We note, however, 
that if part of energy deposited to the medium by the 
traversing parton produces hadrons that remain within 
the angular selection of our jet-peaks, it will not be 
observed in the relative imbalance variable. Thus it is 
possible, that we only measure part of the lost energy 
that got redistributed to hadrons that are softer than 
those considered in this study, or that were radiated 
at larger angles than those included in our jet-cone 
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TABLE I: Jet energy estimate for the near- and away-sides of back-to-back triggered correlations from 200 GeV d+Au data. 
All units are in GeV or GeV/c. For both primary trigger selections listed, the secondary trigger is required to have transverse 
momentum pr G [4, 10] GeV/c. 



Trigger 1 


{Et) T1 


<Pt) T2 


Associated pr 


assoc. near 


assoc,away 

L^iPt 


A(EE T ) 


8 < E T < 10 


8.89 


5.56 


[1.0,10.0] 
[1.5,10.0] 
[2.0,10.0] 
[2.5,10.0] 


4.09±0.06 
3.49±0.06 
2.95±0.06 
2.47±0.06 


4.68±0.06 
3.94±0.06 
3.26±0.06 
2.62±0.05 


2.74±0.09 
2.87±0.08 
3.02±0.08 
3.18±0.08 


10 < E T < 15 


11.74 


5.86 


[1.0,10.0] 
[1.5,10.0] 
[2.0,10.0] 
[2.5,10.0] 


4.14±0.07 
3.59±0.07 
3.09±0.07 
2.65±0.07 


5.75±0.08 
4.88±0.08 
4.02±0.07 
3.32±0.07 


4.26±0.11 
4.58±0.10 
4.95±0.10 
5.20±0.10 



TABLE II: Jet energy estimate for the near- and away-sides of back-to-back triggered correlations from 200 GeV central Au+Au 
collisions. All units are in GeV or GeV/c. For both primary trigger selections listed, the secondary trigger is required to have 
transverse momentum pr £ [4, 10] GeV/c. 



Trigger 1 




(Pt) T2 


Associated pr 


assoc. near 

I^Pt 


assoc,away 

I^Pt 


A(E£ T ) 


8 < E T < 10 


8.87 


5.27 


[1.0,10.0] 
[1.5,10.0] 
[2.0,10.0] 
[2.5,10.0] 


4.68±0.27 
3.81±0.20 
3.16±0.14 
2.61±0.11 


4.82±0.25 
3.86±0.17 
2.94±0.12 
2.13±0.09 


3.45±0.37 
3.55±0.26 
3.82±0.19 
4.08±0.15 


10 < E T < 15 


11.76 


5.70 


[1.0,10.0] 
[1.5,10.0] 
[2.0,10.0] 
[2.5,10.0] 


4.12±0.29 
3.55±0.21 
3.19±0.14 
2.78±0.11 


4.96±0.24 
4.00±0.17 
3.10±0.12 
2.27±0.10 


5.21±0.38 
5.61±0.26 
6.15±0.19 
6.56±0.15 



selection. The dependence of observed energy imbalance 
on the associated hadron pt seems to confirm this idea. 
The relative imbalance is found largest for the highest 
associated hadron px bin studied, reaching the value 
~1.5 GeV/c, and decreasing for the softer associated 
hadrons to approximately ~1 GeV/c. This trend can 
indicate that the jet fragmentation is shifted to the 
lower pt regime (softening), or that the lower energy 
fragmentation products are at larger angles from the 
jet axis (broadening). The later is less likely due to 
consistent jet-like peak widths observed in Au+Au and 
d+Au events. In cither of the two cases, in the limit 
of associated hadron px approaching zero, the value of 
relative A(££t) Au+Au - A(E£: T ) d + Au consistent with 
zero would indicate that full jet energy is recovered. 
Extrapolating our data to very low momenta, albeit 
large uncertainties, seems consistent with this scenario. 

Additional information could be obtained from the rel- 
ative di-jct production rates. For the surface-emission 
scenario the jets and di-jets production rates are deter- 
mined by the surface volume of the fireball. In simplest 
implementation of this scenario the medium in heavy-ion 
collisions is modeled consisting of a completely opaque 
core (leading to full jet absorption) surrounded by a per- 
meable corona (with no jet- medium interactions). In 
such a model a parton (jet) will survive with no energy 
loss or modifications if the underlying hard-scattering 



happened in the corona region and the scattered parton 
did not pass through the core. A Monte Carlo Glauber 
model based on such a model was compared to the early 
di-jet production rates in (Tlj . The size of the permeable 
corona was first determined from tuning to the measured 
single charged hadron Raa results from RHIC. It was 
then used to predict the di-jet survival probability for 
the back-to-back triggers. To minimize systematic un- 
certainties related to the determination of the number 
of binary collisions (N co \\) and the number of participat- 
ing nucleons (-/V part ) in the Glauber model, the double 
ratios measuring conditional di-jet survival probability 

i?, di-jct triggers 

were used. Such a double ratio, -f d +A u " = p A s4i 8 io triggers , 

"AA 

constructed from nuclear modification factors (i?AA, the 
binary scaled ratios of the observables in Au+Au and 
d+Au collisions), reflects any changes in probability to 
find an away-side trigger for each primary trigger ob- 
served in Au+Au data relative to d+Au. In the earlier 
work the value I^+ u Au = 0.20 ± 0.05 for the symmetric 
trigger pairs was found qualitatively consistent with the 
estimates of a simplistic core/corona model [ll|. The di- 
jct survival probability was calculated the same way in 
this work for the new asymmetric trigger data. These 
additional new points for the 0-20% central Au+Au 
data are included in Fig. [11] which reproduces the plot 
from the previous work. We find the conditional sur- 
vival rate ^d+Au " f° r tne asymmetric trigger pairs to be 
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0.26 ± 0.05 for E^ 1 G [8, 10] GeV, and 0.32 ± 0.06 for 
E^ 1 e [10, 15] GcV. We note that the dominant source 
of systematic uncertainties comes from Nu n calculation 
in the Glauber model, which is then fully correlated be- 
tween the data points of the same centrality bin (new val- 
ues) and largely correlated between the values for central 
collisions between the data from different runs. So the 
slight increase in di-jet conditional survival with increas- 
ing trigger pair asymmetry points to higher in-mcdium 
path-lengths, which is consistent with our conclusions 
based on the relative energy imbalance. In other words, 
the new values of I^+Au " indicate that jets triggered by 
those more asymmetric back-to-back pairs are emerging 
from deeper within the "core" of the medium compared 
to a more symmetric pair selection. 




100 200 300 M 

'^part 

FIG. 11: (Color online) Conditional di-jet survival probabil- 
ity measured through di-jet relative suppression rates I^j^ u 
as a function of centrality. The circles show data from [111 ]. 
Filled symbols show new measured ratios for the asymmet- 
ric trigger pairs with primary triggers of [8,10] GeV and 
[10,15] GeV. The band shows the expectation for di-jet surface 
emission ("corona" only) rates from [ill ]. 



VI. SUMMARY 

Jet-medium interactions were studied via a "2+1" 
multi-hadron correlation technique. For both symmetric 
and asymmetric trigger pairs, the distributions of asso- 
ciated hadrons and their spectra show no strong shape 



modifications from d+Au to central Au+Au collisions on 
both near- and away- trigger sides. This is in contrast to 
the di-hadron correlation results with respect to a single 
high-j»T trigger. In addition, no evidence of the near-side 
"ridge" 3 or the away-side "shoulders" Q was observed 
in these "2+1" correlations. The relative total transverse 
momentum imbalance was measured as an excess in the 
difference between the sum of the momentum (or energy) 
for hadrons attributed to the near-side and away-side jet- 
like peaks in Au+Au with respect to the reference d+Au 
data. The relative imbalance of Au+Au over d+Au has 
been shown to be non-zero for the asymmetric trigger 
pair selections, and increasing with the asymmetry of the 
trigger pair. This is consistent with expected medium ef- 
fects for partons probing deeper within the medium. This 
challenges the simplistic implementation of the full ab- 
sorption "core/corona" model, which captured well the 
centrality trend of di-jet production rates for the sym- 
metric triggers. The relative imbalance is found largest 
for the highest associated hadron px bin studied, reach- 
ing the value ~1.5 GcV/c, and decreasing for the softer 
associated hadrons. This trend indicates that the energy 
missing from the away-side peak at higher associated mo- 
menta is converted into softer hadrons. The measured 
relative imbalance for all bins is less than the theoreti- 
cal predictions for total in-medium energy deposition for 
such trigger pairs of ~3 GeV/c in the path-length de- 
pendent energy loss model We note, however, that if 
part of energy deposited to the medium by the traversing 
parton produces hadrons that remain within the angular 
selection of our jet-peaks, it will not be observed in the 
relative imbalance variable. 
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